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Abstract 

Magnetization process of the Gutzwiller wave function is studied accurately 
by a variational Monte Carlo method. We apply it to the one-dimensional 
(ID) and 2D Hubbard models (HM), and to the ID periodic Anderson model 
(PAM) without orbital degeneracy. For the HM, magnetization varies dis- 
continuously to the full moment, as the magnetic field increases. For the 
PAM, the paramagnetic state is unstable against ferromagnetism, although 
the energy reduction thereof is small. 
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Recent experiments under strong magnetic field have renewed the inter- 
ests for heavy-electron systems, such as metamagnetic transitions accompa- 
nied by variations of the Fermi surface and the gap vanishing in the Kondo 
insulators. In this paper, with these in mind, we study properties of the 
Gutzwiller wave function (GWF) |Q] under magnetic field, in particular mag- 
netization process, applying it to the Hubbard model (HM) and the periodic 
Anderson model (PAM). 

The GWF with magnetization m is written as 

|*f(m))=7^(0l*FH), (1) 

where |$p(m)) is a spin polarized one-body function with m = (iVf— N\ )/N a 
(N a and iV a being the numbers of spin a and site, respectively), and V (£) = 
flj 1 — (1 — O n it n i\ ' ^ being the orbital index, and £ a variational param- 
eter which adjusts double occupancy of electrons. Although the GWF is 
a primary many-body wave function, most variational works so far have 
resorted to the so-called Gutzwiller approximation (GA)|l], |2|, ||, |j] in es- 
timating variational expectation values. Here we use a variational Monte 
Carlo (VMC) method for exact evaluation, which tells us how to improve 
the trial wave function. 

First, we discuss the single-band HM, 

n = -t J2 ( c L c j<? + 4*°*) + U J2 ^ n u- ( 2 ) 

<i,j>cr i 

We concentrate on the half- filled band: n = (iVf + N^)/N & = 1, where the 
magnetic-field effect is the most dominant. The energy due to the applied 
field H is given by a Zeeman term: 7"4xt = —g^BHJ2i^i- For this model, 
<3?f("i<) in Eq. (1) becomes a simple Fermi sea, \&F(m)} = n<Tllfc<fc' T c lo-|0)i 
and the effect of magnetic field appears only in magnetization. To obtain 
a magnetization curve, first, we minimize the energy expectation values for 
zero field, EQ(m), with respect to the variational parameters (£ in this case) 
for each m by VMC calculations. Then, we find the relation between m and 
H* by determining m which minimizes Et Q t(H* , m) = -E'o(m) — H*m for a 
given value of H* = gfiBH/2. 

Figure 1 shows the magnetization curves in one dimension (ID) thus 
obtained for some values of U/t with those of two other methods. Under a 
weak field, these three results agree well. However, at some critical field H*, 
m of the GWF (and the GA) saturates discontinuously in contrast to the 
exact one. The jump of m becomes large as U jt increases. The origin of this 
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discontinuity is negative d 2 EQ(m)/dm 2 for m ~ 1, which makes E^ot have 
double minima for H* ~ il*. A similar feature can be seen for large values 
of U/t ( ^ 5) in 2D, as shown in Fig. 2. For small U/t, however, m increases 
smoothly. This behavior of m resembles that of the GA[^] except for the 
existence of the metal-insulator (Brinkman-Rice) transition at U = U c . If 
one uses half ellipse as the density of state, the discontinuity occurs for 
U/U c > 0.44.§ 

In addition to the ID exact result, || quantum Monte Carlo calcula- 
tions in infinite dimension^ concluded no discontinuity in the magnetiza- 
tion curve. We consider that this discrepancy is attributed to the onsite 
nature of the correlation factor in the GWF; an intersite correlation fac- 
tor, spin dependent in this case, will encourage the electron transfer, which 
becomes more important as m increases. 

Incidentally, for the ID t-J model the GWF does not have a jump in 
the magnetization curve for any value of J/t and well reproduces the exact 
magnetic properties for the half filling or the supersymmetric case.Q] This 
is because the intersite correlation effect is introduced through the canonical 
transformation from the HM to the t-J model. || 

Having restricted n to the half filling, we have confirmed that the ten- 
dency mentioned here holds also for n < 1. 

Next, we consider the ID periodic Anderson model without f-level de- 
generacy, 



ka ka 

+E £ f/L/^ + C7 E^X'i- (3) 

ka j 

For simplicity, we put = —2tcosk, Ef and Vk{= V) constant, and fix 
two parameters at typical and convenient values: V/t = 0.5 and U/t = oo. 
Furthermore, as for the Zeeman term we assume g c = 9i- 

For this Hamiltonian, we consider the Gutzwiller-projected hybridized 
band state: Eq. (1) with P f (0) and 

IMm)> = II II H^4. + sm^/lj |0), (4) 

a k<k° 

where we choose as a spin-independent noninteracting form: 

2V 

tan^fc = -j^^^^=^^^^= (5) 

et-e k ± J (ii - e k f + AV 2 



3 



for the upper (— ) and lower (+) bands, respectively, if and V are variational 
parameters, which adjust the band form. This wave function is an extention 
of Yosida's singlet state for the Kondo problem, Q and has been extensively 
studied for the cases without magnetic field. ||i~Cfl 

In Fig. 3 the energy expectation values without field are depicted. We 
have calculated for metallic (n = 1.6) as well as insulating (n = 2.0) 
cases. In each case, EQ{m)/t varies slightly for small m and then in- 
creases abruptly with increasing m. This aspect is similar to the exact- 
diagonalization results. [|ll]] Meanwhile, the optimized state has finite spin 
polarization, as well as dEo(m) / dm\ m= o > and d 2 Eo(m)/dm 2 \ m= o < 0. 
Also for a finite U (U, V, Ef = l.Ot, 0.5t, — l.Ot), we have confirmed that the 
above derivatives have the same signs with those for infinite U ; namely the 
GWF is unstable against a magnetic order. These aspects agree with the 
Gutzwiller approximations, |3], Q but does not with an exact analysis for the 



symmetric case, [12] and a study in infinite dimension, [13] where the ground 
state is paramagnetic. ^From these facts we observe that the ground state of 
the ID nondegenerate PAM, either paramagnetic or partially ferromagnetic, 
is rather fragile; actually the energy differences in Fig. 3 are the order of 
O.Olt. Stable paramagnetism may need f-level degeneracy. Q 

Magnetization curves obtained from Fig. 3 are shown in Fig. 4. There 
is no sign of metamagnetism. It also remains a future problem how the gap 
for the insulating case changes as a function of various parameters. 

In conclusion, although the VMC results support the GA qualitatively, 
improvements of the trial function are needed for advanced discussions. To 
this end, introduction of a spin-dependent one-body part as well as of an 
intersite (RKKY-like) or off-diagonal correlation factors is important. 
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Figure captions 

Fig. 1 Magnetization as a function of magnetic field for the ID Hubbard 
model. The results of the GWF are shown with those of the GA and the 
Bethe Ansatz (Exact) f§ for U/t = 3. The system used for the GWF has 90 
sites with closed-shell condition. For the GA the ID cosine band is assumed. 

Fig. 2 Magnetization as a function of magnetic field of the Hubbard model 
for the 2D square lattice by the GWF. We use systems of L x L (L = 
10, 12, 14) sites with the antiperiodic-antiperiodic boundary condition and 
closed shell condition. 

Fig. 3 Energy expectation values without magnetic field as a function of 
magnetization for ID periodic Anderson model. Minimal values for Ef/t = 
—1.0 are plotted by dotted lineas as a guide. The system of N a = 50 is used 
with closed shell condition. 

Fig. 4 Magnetization as a function of magnetic field of the ID periodic 
Anderson model for some parameter values. For comparison, the data for 
U/t = 0, V/t = 0.5 and -Ef/t = —1.0 are plotted by dash-dotted lines. 
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